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Thermal fluctuations in time-dependent quantum processes are treated by a constant-temperature 
generalization of Wigner's formulation of quantum mechanics in phase space. To this end, quantum 
Nose-Hoover dynamics is defined by generalizing the Moyal bracket. Computational applications of 
the formalism, together with further theoretical developments, are discussed. 



Many (equivalent) ways to represent quantum mechan- 
ics in phase space [l| have been developed since the orig- 
inal paper of Wigner Although the phase space is 
comprised of the classical coordinates of the system (r, p) , 
positions and momenta respectively, it is nevertheless 
quantum since the laws of motion and the statistical con- 
straints of the theory are designed so that standard quan- 
tum averages can be reproduced. From a formal point of 
view, we may note that, within all the different formu- 
lations of quantum mechanics in phase space, Wigner's 
provides perhaps the most simple one. Wigner's formula- 
tion of quantum mechanics has been recently used to in- 
vestigate quantum states of the electromagnetic field Q , 
quantum-to-classical correspondence both in the elec- 
tronic transport properties in nanowires [4[ and in driven 
nonlinear nanomechanical resonators [5] . 

Wigner's approach to quantum mechanics can be for- 
mulated by means of the Moyal bracket jy] . This provides 
an algebraic route to the formulation of quantum me- 
chanics in phase space d, 0] which has also been exploited 
in what seems to be an emerging branch of quantum the- 
ory, known as deformation quantization [8(. As a mat- 
ter of fact, Moyal's bracket is obtained within the defor- 
mation quantization approach upon deforming the well- 
known classical Poisson structure in order to introduce 
a generalized bracket in phase space expressing quan- 
tum laws. We note that a strategy similar to that used 
in deformation quantization, i.e., the generalization of 
brackets, has been adopted very recently to introduce a 
formalism to derive generalized equations of motion and 
analyze the statistical mechanics of non-Hamiltonian sys- 
tems both in the classical @ and in the quantum con- 
text fToL rTTj| . In particular, it has been shown [13] how to 
formulate constant-temperature dynamics for quantum- 
classical systems so that a reduced number of classical 
bath degrees of freedom can be used to simulate relax- 
ation dynamics. This required to combine the method 
of constant-temperature dynamics [l3j (known as Nose- 



Hoover dynamics) used in computer simulations of clas- 
sical systems with the quantum-classical dynamics [ill ]. 
In classical molecular dynamics the Nose-Hoover equa- 
tions [Hi], and their variants [3], sample the canonical 
ensemble very efficiently and they can all be formulated 
as non-Hamiltonian phase space flows. 

In many instances, ranging from quantum phase tran- 
sitions [HI to lattice gauge theory simulations 16], one is 
interested in calculating quantum averages in the canon- 
ical ensemble. Hence, one could enquire whether, within 
the Wigner formulation of quantum mechanics in phase 
space, a non-Hamiltonian generalization of the Moyal 
bracket can be introduced and used to investigate the ef- 
fect of thermal fluctuations in the applications mentioned 
above [3, \4, ki- 
ln this letter we introduce constant-temperature Nose- 
Hoover dynamics [H| in quantum phase space. To this 
end, we generalize the Moyal bracket [H, HI in a non- 
Hamiltonian fashion [|| [l(| [Hj]. Such a generalization 
allows one to devise systematic approaches to the cal- 
culation of quantum corrections to equilibrium or time- 
dependent properties taking explicitly into account ther- 
mal fluctuations. 

Given the von Neumann density matrix p, the quan- 
tum average of an arbitrary operator x can be calculated 
in phase space by introducing the Wigner transform of 
both the operator and the density matrix: Xw{ r iP) = 
W(x), fw(r,p) = {2nh)- N W{p), where W denotes the 
operator which realizes the Wigner transform [2;], N is 
the number of degrees of freedom of the system, and (r, p) 
are the phase space coordinates. The phase space quan- 
tity fw{ r ,p) is known as the Wigner function. Then, the 
quantum average (x(i)) = Tr(x(t)p) can be calculated as 



(X) 



d rd p Xw(r,P,t)fw(r,p) ■ 



(1) 



The law of motion for an arbitrary operator x can be 
found by Wigner-transforming the Heisenberg equation 
of motion: 



9tXw(r,P,t) = t {H w ,xw{t)} 



M ' 



(2) 
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where d t = d/dt, Hy/ = W(H) is the Wigner equivalent 
of the Hamiltonian operator H. The right hand side of 
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Eq. Q defines the Moyal bracket, which for two arbitrary 
(quantum) phase space functions, Xw ano - Xw> can ^ e 
written as 



{Xw,Xw} M = X^exp 



Xw ex P 



(3) 



where we have indicated the phase space point through 
the compact notation x — (£1,2:2) = {r,p), and di stands 
for d/dxi. Note that we are using Einstein's convention 
of summation over repeated indices, i,j which go from 1 
to 2N, and we have introduced the antisymmetric matrix 



1 
1 



(4) 



Equation ([3]), which recasts the Moyal bracket in matrix 
form, is central to the formalism that will be presented 
in the following. Defining the operator 



M = H\ 



tr 



d 3 _ p -fd i Bt 3 d ] 



(5) 



the propagator and the time evolution of observables can 
be written as 



Xw(r,p,t) = exp 



it 



M 



Xw(r,p) ■ 



(0) 



Together Equations {1} and ([6]) specify how to calculate 
quantum averages in phase space. By means of repeated 
partial integrations, the time dependence in {!]) can be 
transferred from the observable xw to the Wigner func- 
tion /w, eventually recovering Wigner's equation of mo- 
tion for the latter |2j. For standard quantum dynamics, 
this amounts to substitute fw for xw in Eq. |6|) and hav- 
ing a minus sign in the propagator (exp —j^M ). In a 

nutshell, this expresses the essence of the Wigner-Moyal 
formulation of quantum mechanics. 

The matrix form of the Moyal bracket, introduced in 
Eq. (|3|), provides a mathematical structure that can be 
naturally generalized for defining a Nose-Moyal bracket 
with a non-Hamiltonian character. In order to introduce 
such a Nose-Moyal bracket, we extend the quantum phase 
space by adding two Nose variables r),p v , with fictitious 
mass m„, as it is done in the classical Nose- Hoover molec- 
ular dynamics approach. Hence, the phase space point 
coordinates are given by x — (r,rj,p,p n ). Correspond- 
ingly, we introduce a quantum extended Nose Hamilto- 



V 2 

H w {x) = H w (r,p) + -p- + gk B Tr, 



2m, 



(7) 



where g is the number of degrees of freedom, kg is Boltz- 
mann's constant, and T is the fixed temperature of the 
canonical ensemble. In the extended phase space, we can 



now generalize the symplectic matrix B c to the antisym- 
metric tensor [t| ITU Il2l| 
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(8) 



and re-define the Moyal bracket in a way suitable for the 
derivation of Nose-Hoover dynamics. In analogy to the 
definition (|5|), we can use the antisymmetric tensor S N 
to introduce an operator in the extended quantum phase 
space 



M 



ttN 



e 2 



r di By d j 



(9) 



Equation © represents the fundamental Nose-Moyal op- 
erator for achieving Nose-Hoover dynamics in quantum 
phase space. The time-propagation of the observable can 

— >N 

be re-defined by using the Nose-Moyal operator M. in 
place of Ai in Eq. ^ . Again, by means of repeated par- 
tial integrations, the time dependence can be transferred 
from xw to fw- However, in the extended phase space 
one finds 



fw{x,i) = exp 



it—>N,t 



fw(x) 



(10) 



with 



M = 



rrN 



exp 



exp 



(11) 



This means that, under Nose dynamics, the stationary 
Wigner function, fw,e, obeys the following equation 



iL N f w , t 



« JW,e 



x H, 



E 

n=3,5,7,. 



1 



IW.e 



where 



(fyBg) 3 H W = Np./m, 



(12) 

(13) 
(14) 



are the Nose phase space compressibility and the Nose- 
Liouville operators, respectively. Equations (TiT)]) and (Tl2]) 
contain the full quantum corrections to the dynamics of 
all the variables in the extended phase space. They define 
Nose-Hoover dynamics in a completely quantum fashion. 
One of the main interests of such a generalization lies in 
the fact that, as in the classical case, just a pair of ad- 
ditional variables, namely the Nose coordinates (t?,^), 
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allows to represent thermal fluctuations (and therefore 
the process of relaxation toward thermodynamical equi- 
librium) of a quantum system by introducing a suitable 
non-Hamiltonian dynamics in the extended phase space. 
This is to be confronted with a standard Hamiltonian 
formalism that would require coupling the physical coor- 
dinates, (r,p), to a bath composed by an infinite number 
of degrees of freedom. 

Once a stationary expression of the Wigner function 
has been found, equilibrium quantum averages that in- 
clude thermal fluctuations in the dynamical evolution 
(and not just in the initial conditions) can be calculated 
by propagating observables according to 



dtXw = ~T H wsin 

h 



Xw 



(15) 



Equation (|15[) shows how to propagate a dynamical vari- 
able undergoing Nose- Wigner time evolution in quan- 
tum phase space taking into account quantum effects 
also in the dynamics of the Nose coordinates. In or- 
der to make clear the characteristics of the quantum 
corrections, one can expand in a Taylor series the sine 
in the right hand side of Eq. (fT5| . Upon realizing that 
d 2 H^/dif = d 3 H^/d P 3 = d^H^./dpl = and that all 
mixed derivatives d 2 H^/dxidxj = 0, for i ^ j, Eq. (TH 
can be rewritten as 



dtXw = iL xw 



E 



n=3,5,7,. 



+ H W 




rrN v 

Hw d^dr Xw 



dp v \ di] 1 1 dp 



Xw 



(16) 



We remark that, since thermal effects are described 
within the dynamics, one can also study relaxation pro- 
cesses by using a non-equilibrium Wigner function and 
evolving the dynamical variable according to Eq. (j 1 5jl : 
the dissipation will be forced by the non-Hamiltonian 
coupling to the Nose variables, which mymic the ther- 
modynamic of a thermal bath. 

The analysis of the stationary Nose-Hoover distribu- 
tion in quantum phase space can be simplified by per- 
forming a classical approximation on the dynamics of the 
Nose variables alone and keeping all quantum corrections 
in the dynamics of the coordinates of the physical system 
(r, p). In fact, we remind that, when Nose-Hoover dy- 
namics is implemented within molecular dynamics com- 
puter simulations, a mass ~ N x m is used in or- 
der to achieve a weak coupling to the Nose "bath" . In 
such a way, while the temperature of the system is con- 
trolled, the equilibrium dynamical properties of the phys- 
ical coordinates (r, p) are not significantly modified. In 
the present context, a small parameter p = y^m/m n is 



naturally found in the theory so that one is allowed to 
perform a classical limit on the Nose coordinates {mp^). 
In this way, a quantum-classical description, along the 
lines described in [l7l ]. naturally arises. Therefore, within 
a quantum-classical approach, one could disregard all the 
quantum corrections on the evolution of the Nose vari- 
ables in the left hand side of Eq. (fT2"| . Upon assuming a 
standard form for the Hamiltonian of the physical degrees 
of freedom Hw{r,p) = (p 2 /2m) + V(r), the Nose-Hoover 
equation for the stationary Wigner distribution function 
becomes 



E 







n 




d r ■ ~£) p 


fw,e 



n=3,5,7, 



(17) 



where d r — d/dr and d p = d/d p . Moreover, in order 
to calculate quantum averages of functions of (r, p) alone 
we need to calculate the average of Eq. p^|) over the 
now classical Nose variables. This turns out to be iden- 
tical to the equations first proposed by Wigner 0], who 
showed how to obtain quantum corrections in the canon- 
ical ensemble in terms of an expansion of his distribution 
function in even powers of U [2| : 



E 

n=0 



An) 
JW,e 



(18) 



The order zero solution, fjjp, to Eq. (fT?]) is [|[ 

f$ot5(H%)e-J K " dt . (19) 

When averaging over 77, p n , f$) becomes the Boltzmann's 
weight for the (r, p) variables. Moreover, averages of odd 
functions of the Nose variables over fff are zero. Hence, 
the analysis of the correction terms, under the approxi- 
mations of no quantum effects on the Nose variables, can 
proceed as originally shown by Wigner 0] • The higher or- 
der correction terms will all contain the zero-order term: 



,{n) 
J\V,e 



ft" fwefwe' w ith n > 2. In principle, such cor- 



rection terms could be used in order to calculate aver- 
ages over an ensemble of dynamical trajectories including 

thermal fluctuations by means of Nose-Hoover dynam- 

— >n 

ics as defined through the operator M. introduced in 
Eq. ©. 

Previous attempts to formulate quantum Nose-Hoover 
dynamics 18] run into the problem of the nonergodic- 
ity of the dynamics. However, in the present case, ex- 
ploiting the antisymmetric matrix structure of Moyal 
bracket, introduced Eqs. ([5]) and ([9]), one can easily 
formulate more sophisticated thermostatting methods, 
such as Nose- Hoover chain dynamics 14], in order to, at 
least, tame eventual ergodicity problems with dynamics 
in quantum phase space. 

In addition to the Nose-Hoover dynamics in quantum 
phase space, we would also like to discuss some interest- 
ing theoretical developments which arise from the matrix 
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structure of the Moyal bracket, which has been intro- 
duced in Eqs. ([5]) and ©. Such a matrix structure allows 
suggests a connection between the approach of deforma- 
tion quantization and a formalism proposed by Wein- 
berg in order to address non linear effects in quantum 
mechanics Weinberg's formalism has been slightly 

generalized and expressed in matrix form in [Io| . Having 
introduced a generalized matrix structure for the Moyal 
bracket, Weinberg's non-linear formalism can be re-cast 
within a deformation quantization approach: one just 
needs to define more general Moyal brackets by means 
of an antisymmetric tensor 3 = B[fw{%)]> which can 
be a functional of the Wigner distribution function itself. 
To the authors's knowledge such mathematical structures 
have not been been studied so far. In such a way, one 
has at disposal a deformation quantization formalism for 
addressing eventual non-linear phenomena in quantum 
mechanics. 

From a practical point of view, the main justifica- 
tion for the continuous attempts to modify or general- 
ize the formalism of quantum mechanics lies in the fact 
that the devolopment of a general method to simulate 
quantum dynamics still remains an open problem fac- 
ing the physics community. The availability of such a 
method, of course, defines the very possibility of pre- 
dicting and calculating quantities of interest in many 
physico-chemical processes in condensed matter systems 
(including biological ones). Typically, one desires to for- 
mulate quantum properties in a framework as close as 



possible to the classical one so that computational meth- 
ods, which are known to be efficient in classical mechan- 
ics, can be adapted to study quantum time-dependent 
processes. For example, the path integral approach to 
quantum mechanics provides an isomorphism between 
quantum and classical systems which precisely permits to 
devise such computational schemes. However, the above 
mentioned isomorphism cannot be easily exploited, to ad- 
dress quantum-time dependent properties in an efficient 
way. Other routes, such as the one we have sketched 
here, must be searched. As a matter of fact, there are 
various computational schemes that, exploiting the the- 
ory of stochastic processes [2(1, attempt to simulate ei- 
ther full quantum dynamics [2l| or quantum-classical dy- 
namics in phase space (22j . Deformation quantization is 
one of such routes and we believe that, by generalizing 
the Moyal bracket, we have taken a first step toward the 
pursue of a mathematical language for developing novel 
computational approaches. Typically, the introduction 
of the Nose- Moyal bracket in quantum phase space, and 
the corresponding possibility of treating thermal fluctu- 
ations through non-Hamiltonian dynamics, opens a new 
route for treating open quantum systems [201 ] . 
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